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We study the entanglement in the quantum Heisenberg XY model in which the so-called W 
entangled states can be generated for 3 or 4 qubits. By the concept of concurrence, we study the 
entanglement in the time evolution of the XY model. We investigate the thermal entanglement 
in the two-qubit isotropic XY model with a magnetic field and in the anisotropic XY model, and 
find that the thermal entanglement exists for both ferromagnetic and antiferromagnetic cases. Some 
evidences of the quantum phase transition also appear in these simple models. 
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PACS numbers: 03.65.Ud, 03.67.Lx, 75.10.Jm. 
I. INTRODUCTION 

Quantum entanglement has been studied intensely in 
recent years due to its potential applications in quantum 
communication and information processing pi such as 
quantum teleportationB] , superdense coding |1, quan- 
tum key distribution [QfTand telecoloning M. Recently 
Diir et al. H found that truly tripartite pure state entan- 
glement of three qubits is either equivalent to the maxi- 
mally entangled GHZ state M or to the so-called W state 



|W0 = -7=(|100) + |010> + |001)). 



(1) 



For the GHZ state, if one of the three qubits is traced 
out, the remaining state is unentangled, which means 
that this state is fragile under particle losses. Oppositely 
the entanglement of W state is maximally robust under 
disposal of any one of the three qubits H . 

A natural generalization of the W state to N qubits 
and arbitrary phases is 



\W N ) = -We iei |1000...0) + e l92 |0100...0) 
v N 



e i03 10010. ..0) + • • • + e i6N |0000...1)). (2) 

For the above state | Wjy), the concurrences ||D between 
any two qubits are all equal to 2/N and do not depend 
on the phases Bi(i = 1,2...N). This shows that any two 
qubits in the W state are equally entangled. Recently 
Koashi et al. || shows that the maximum degree of en- 
tanglement (measured in the concurrence) between any 
pair of qubits of a iV-qubit symmetric state is 2/N. This 
tight bound is achieved when the qubits are prepared in 
the state \Wn). 

The Heisenberg interaction has been used to imple- 
ment quantum computer fp_ 0B - I t can be realized inquan- 
tum dots |l(J , nuclear spins [|ll[] , electronic spins [G3 and 
optical lattices [|l3|. By suitable coding, the Heisenberg 
interaction alone can be used for quantum computation 



Here we consider the quantum Heisenberg XY model, 
which was intensively investigated in 1960 by Lieb, 
Schultz, and Mattis fllq] . Recently Imamoglu et al have 
studied the quantum information processing using quan- 
tum dot spins and cativity QED [fl6f and obtained an 
effective interaction Hamiltonian between two quantum 
dots, which is just the XY Hamiltonian. The effective 
Hamiltonian can be used to construct the controlled- 
NOT gate fuj. The XY model is also realized in the 
quantum-Hall system |l7j and in cavity QED system |i"§| ] 
for a quantum computer. 

The XY Hamiltonian is given by [JL5| 
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where S a — a a /2 (a — x, y, z) are spin 1/2 operators, a a 
are Pauli operators, and J > is the antiferromagnetic 
exchange interaction between spins. We adopt the peri- 



odic boundary condition, i.e., S 
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One role of the XY model in quantum computation 
is that it can be used to construct the swap gate. The 
evolution operator of the corresponding two-qubit XY 
model is given by 



U(t) = exp [-iJi(af of + crla^/2} . 
Choosing Jt = ir/2, we have 



U(£) |00) = |00>^(^)|11> 



(4) 



111) 



u (I/) |01> = _l|10> ' u (h) |10> = ^ |01> - (5) 



The above equation shows that the operator U (jj) acts 
as a swap gate up to a phase. Another gate ^/swap which 
is universal can also be constructed simply as U {jj) ■ A 
swap gate can be realized by successive three C-NOT 
gates |EJ , while here we only need one-time evolution of 
the XY model. This shows that the XY model has some 
potential applications in quantum computation. 

The entanglement in the ground state of the Heisen- 
berg model has been discussed by O'Connor and Woot- 
ters |20|. Here we study the entanglement in the XY 
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model. We first consider the generation of W states in 
the XY model. It is found that for 3 and 4 qubits, the W 
states can be generated at certain times. By the concept 
of concurrence, we study the entanglement properties in 
the time evolution of the XY model. Finally we discuss 
the thermal entanglement in the two-qubit XY model 
with a magnetic field and in the anisotropic XY model. 



II. SOLUTION OF THE XY MODEL 

With the help of raising and lowering operators 
S*±iS%, the Hamiltonian H is rewritten as ( J = 1) 



1 N 



(6) 



Obviously the states with all spins down |0)® or all 
spins up 11)®^ are eigenstates with zero eigenvalues. 

The eigenvalue problem of the XY model can be ex- 
actly solved by the Jordan- Wigner transformation pl| . 
Here we are only interested in the time evolution prob- 
lem and in the 'one particle' states (N — 1 spins down, 
one spin up), 
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The eigenequation is given by 
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(7) 



(8) 



n=l 
N 



= E k J2a k , n (T+\0)^ N . 

n=l 

Then the coefficients ak, n satisfy 
The solution of the above equation is 



(9) 



a k , n = exp 
E k = cos 



i2-Knk\ , 

—) (k=l...N), 

2nk 
If 



(10) 
(11) 



where we have used the periodic boundary condition. 
So the eigenvectors are given by 



N 



v n—1 ^ ' 



(12) 



which satisfy (k\k') — <5fcfc'.It is interesting to see that all 
the eigenstates are generalized W states (Eq.(||)). 



Note that the XY Hamiltonian H commutes with the 
operator 



= af N =a x ®a x 



then the state 



n—1 



(13) 



(14) 



are also the eigenstates of H with eigenvalues 
cos (2nk/N) . 

Now we choose the initial state of the system as 
a i |Q)® , and in terms of the eigenstates \k), it can be 
expressed as 



1 N 
|*(0)> = -=V exp 



-i2nk 
N 



\k). 



(15) 



The state vector at time t is easily obtained as 



\*(t))=J2.bn(t)<7+\0)® N , 



where 



1 N 



)k/N -itcos(2-nk/N). 



(16) 



(17) 



fc=i 



If we choose the initial state as a 1 \1)® N , then the wave 
vector at time t will be J2n=i & n(*)°n l 1 )®^- 

III. GENERATION OF W STATES 

From Eq.(|I7|), the probabilities at time t for state 
a+\0)® N is obtained as 



P(n,N,t) = \b n {t)\ 2 



(18) 



For N = 2, it is easy to see that the probability 
P(l,2,i) = cos 2 i , P{2,2,t) = sin 2 t. The state vector 
at time t is 



\^{t)) = cost|10) - isint|01). 



(19) 



When t — 7r/4, the above state is the maximally entan- 
gled state. 

Now we consider the case N — 3. The probabilities are 
analytically obtained as 



P(l,3,t) 



5 + 4 cos ( -t 



P(2,3,i) = P(3,3,t) = - 



2 - 2 cos I -t 



(20) 
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Fig. 1(a) gives a plot of the probabilities versus time. It is 
clear that there exist some cross points of the probabili- 
ties. At these special times the probabilities P(n, 3, t) are 
all equal to 1/3, which indicates the W states are gener- 
ated. From Eq. (|2C|) , we see that if the time t satisfies the 
equation 



'3 \ 1 



(21) 



the probabilities are same. The solution of Eq. (|2l|) is 



t = 



~9~ 
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(" = 0,1,2,. 



(22) 



Explicitly at these time points, the corresponding state 
vectors are 



As seen from Fig. 2(a), there also exists some cross 
points, which indicates the 4-qubit W states are gen- 
erated. The probabilities are same when 



t n = - + 2mr, 
t' n = y + 2mr(n = 0,l,2,...) 
Explicitly the 4-qubit W states are 



(25) 



!*(*»)) 
l*CO) 



- fllOOO) - i|010Q) - |0010) - i|0001» , 

- (|1000) + i|0100) - |0010> + i|0001» . 

(26) 



|*(*n)> = ^(|100)+e^|010) + e^|001)), 
|*(f B )) = ^(|100)+e*|010)+e*|001))(23) 
which are the generalized W state for TV = 3. 




FIG. 1. Time evolution of the probabilities and concur- 
rences for iV = 3. (a)The probablity plus 1 for n — 1 (solid 
line) and n = 2, n = 3 (dotted line); (b)The concurrences 
Ci 2 (t),Ci 3 (t) (solid line) and C 23 (f) (dotted line). 




FIG. 2. Time evolution of the probabilities and concur- 
rences for N = 4. (a)The probablity plus 1 for n = 1 (solid 
line), n = 3 (dotted line) and n = 2, n — 4 (dashed line); 
(b)The concurrences Cvi{£) (solid line) and C23,(t) (dotted 
line) . 

Can we generate W states for more than 4 qubits in 
the XY model? Fig. 3 (a) shows that there is no cross 
points for N = 5 . Further numerical calculations for 
long time and large TV show no evidence that there exist 
some times at which the W states can be generated. 



For the case N = 4, the probabilities are given by 



P(l,4,t) = cos 4 (£) , P(3,4,i)=sin 4 Q 
P(2,4,t) = P(4,4,i) = isin 2 i, 



(24) 
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defined as 



C12 = max {Ai - A 2 - A 3 - A 4 , 0} 



(27) 
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FIG. 3. Time evolution of the probabilities and concur- 
rences for N = 5. (a)The probablity plus 1 for n — 1 
(solid line), n = 3, n = 4 (dotted line) and n — 2, n = 5 
(dashed line) ;(b) The concurrences Ci2(i) (solid line) and 
C 3 a(t) (dotted line). 

We see that the W states appear periodically for 3 
and 4 qubits. In order that a certain state occurs period- 
ically in a system, a necessary condition is that the ratio 
of any two frequencies available in the system is a ratio- 
nal number. From Eq.([fl]) it is easy to check that the 
necessary condition is satisfied for 2, 3, 4, and 6 qubits. 
For 6 qubits the corresponding probabilities evolve pe- 
riodically with time. The numerical calculations show 
that there exists no cross points, i.e., we can not creat 
6-qubit W state. However some states close to the W 
state will occur repeatedly and these states may be used 
for quantum computation. 

The 3-qubit and 4-qubit W states are readily generated 
by only one-time evolution of the XY system. This idea 
is similar to the concurrent quantum computation p2| ] in 
which some functions of computation are realized by only 
one-time evolution of multi-qubit interaction systems. 

The W entangled states can be generated by other 
methods, such as coupling N spins with a quantized elec- 
tromagnetic field. However here we only use the interac- 
tion of N spins themselves and do not need to introduce 
additional degree of freedoms. 



IV. TIME EVOLUTION OF ENTANGLEMENT 

We first briefly review the definition of concurrence 
M. Let pi2 be the density matrix of a pair of qubits 1 
and 2. The density matrix can be either pure or mixed. 
The concurrence corresponding to the density matrix is 



where the quantities Ai > A 2 > A3 > A4 are the square 
roots of the eigenvalues of the operator 



(28) 



The nonzero concurrence implies that the qubits 1 and 
2 are entangled. The concurrence Cyi = corresponds 
to an unentangled state and C12 = 1 corresponds to a 
maximally entangled state. 

We consider the entanglement in the state |^(£))(|l6j). 
By direct calculations, the concurrence between any two 
qubits i and j are simply obtained as 



C ij (t)=2\b i (t)b j (t)\. 



(29) 



The numerical results for the concurrence are shown in 
Fig.l(b), Fig.2(b) and Fig.3(b). 

For TV = 3, Fig. 1(b) shows that the entanglement is pe- 
riodic with period 47r/3. At times 4n7r/3(n = 1,2,3, ...), 
the state vectors are disentangled and become the state 
1 100) up to a phase. The concurrences of C 12(f) an d 
Cis(t) are same, and have two maximum points in one 
period, while the concurrence C\^(t) has only one maxi- 
mum point. Fig. 2(b) shows the concurrences for N = 4. 
They are periodic with period 27r. In one period there are 
two unentanglemcnt points, t — it, 2tt. For both concur- 
rences C\2(t) and C*23(t), there are two maximum points 
in one period. If we choose large N (see Fig. 3(b) for 
N = 5), there exists no exact periodicity for the entan- 
glements of two qubits. From the time evolution of the 
concurrences we can see clearly when the system becomes 
disentangled and when the system maximally entangled. 



V. THERMAL ENTANGLEMENT 

Recently the concept of thermal entanglement was in- 
troduced and studied within onc-dimensional isotropic 
Heisenberg model (^3|. Here we study this kind of en- 
tanglement within both the isotropic XY model with a 
magnetic field and the anisotropic XY model. 



A. Isotropic XY model with a magnetic field 

We consider the two-qubit isotropic antiferromagnetic 
XY model in a constant external magnetic field B, 



B 



H=-K+ erf) + J (<J+<J 2 + a+a l ) 



(30) 



The eigenvalues and eigenvectors of H are easily obtained 
as 



£T|00) = -fl|00),£r|ll) = B|ll) 



m^) = ±.m ± ) 



(31) 



4 



TTadOl) ± 1 10) ) are maximally entangled 



where |\I' ± ) 
states. 

The state of the system at thermal equilibrium is 
p(T) = exp(-jg,) /Z, where Z =Tr[exp(-^)] is the 
partition function and k is the Boltzmann's constant. As 
p(T) represents a thermal state, the entanglement in the 
state is called thermal entanglement |^3| . 

In the standard basis, {|00), |01), 1 10) , |11)}, the den- 
sity matrix p(T) is written as {k = 1) 



P(T) 



1 



2(cosh ^ + cosh ^) 



(32) 
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Vo 





cosh ^ 
— sinh y 






— sinh ^ 
cosh y 









S- I 

eT / 



From Eqs.(27),(|28|) and (32), the concurrence is given by 



C = max 



sinh ^ — 1 
cosh y + cosh y 



Then we know C = if sinh -yr < l,i.c. 



(33) 



there is a 



critical temperature 



J 



arcsinh(l) 



1.1346J, 



(34) 



the entanglement vanishes for T > T c . It is interesting to 
see that the critical temperature is independent on the 
magnetic field B. 

For B = 0, the maximally entangled state is 
the ground state with eigenvalue — J. Then the maxi- 
mum entanglement is at T — 0, i.e., C — l.As T in- 
creases, the concurrence decreases as seen from Fig. 4 due 
to the mixing of other states with the maximally entan- 
gled state. For a high value of B (say B — 1.2), the state 
|00) becomes the ground state, which means there is no 
entanglement at T = 0. However by increasing T, the 
maximally entangled states |\E' ± ) will mix with the state 
|00), which makes the entanglement increase (see Fig. 4). 
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FIG. 4. The concurrences versus temperature for different 
magnetic fields. The parameter J is set to one. 

From Fig. 5 we see that there is a evidence of phase 
transition for small temperature by increasing magnetic 
field. Now we do the limit T — > on the concurrence 
(p3|), we obtain 



lim C 

lim C 
lim C 



= 1 for B < J, 

= - for B = J, 
2 

= for B > J. 



(35) 



So we can see that at T = 0, the entanglement vanishes 
as B crosses the critical value J. This is easily understand 
since we see that if B > J, the ground state will be the 
unentangled state 1 00) . This special point T = 0, B = J, 
at which entanglement becomes a nonanalytic function 
of B 1 is the point of quantum phase transition |24||. 
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FIG. 5. The concurrences versus magnetic field B for dif- 
ferent temperatures. The parameter J is set to one. 

It should be pointed out that the results of thermal en- 
tanglement in the present isotropic XY model is qualita- 
tively the same as but quantitatively different from that 
in the isotropic Heisenberg model (23). An important 
conclusion is that the concurrences are the same for both 
positive J and negative J in the XY model. That is to 
say, the entanglement exists for both antiferromagnetic 
and ferromagnetic cases. In contrary to this, for the case 
of two-qubit Heisenberg model, no thermal entanglement 
exists for the ferromagnetic case. 



B. Anisotropic XY model 

Now we consider the two-qubit anisotropic antiferro- 
magnetic XY model which is described by the Hamilto- 
nian jl5j 



J 



[(1+7K^ + (1-7K^] 



(36) 



where 7 is the anisotropic parameter. Obviously the 
eigenvalues and eigenvectors of the Hamiltonian H a is 
given by = ±J\V ± ) and = ±J7|<& ± ), 

where |«& ) = ^(|00) ± |11)). Then the four maximally 
entangled Bell states are the eigenstates of the Hamil- 
tonian H a . Although the anisotropic parameter can be 
arbitrary, we restrict ourselves on < 7 < 1. The param- 
eter 7 = and 1 correspond to the isotropic XY model 
and Ising model respectively. Thus the anisotropic XY 
model can be considered as a interpolating Hamiltonian 
between the isotropic XY model and the Ising model. 
The anisotropic parameter 7 controls the interpolation. 



by 



The density matrix p(T) in the standard basis is given 



P(T) 



( cosh ^ 



V-sinh^ 



2(cosh J +cosh i?) 

-sinh^N 

- sinh 







coshy :„,,,> — 

— sinh y cosh y 







cosh ^ / 



(37) 



The square root of the eigenvalues of the operator Q12 
are — — , e f ' — — ttt an d —, — , e , 7 — : j 7 , ■ Then from 

2(c osh T77+cosh -Tjr ) 2(cosn ^r+cosh ~^~) 

Eq. (p7[) , the concurrence is given by 



„ / sinh i — cosh ^ 

C = max ^ 

V cosh 4p 



- COSh 



(38) 



As we expected Eq. (g8|) reduces to Eq. (|33|) with B = 
when 7 = 0. When 7 = 1, the concurrence C = 0, 
which indicates that no thermal entanglement appears 
in the two-qubit Ising model. In this anisotropic model, 
the concurrences are the same for both positive J and 
negative J , i.e, the thermal entanglement is the same 
for the antiferromagnetic and ferromagnetic cases. The 
critical temperature T c is determined by the nonlinear 
equation 

sinh — = cosh — - , 
T T ' 

which can be solved numerically. 

In Fig. 6 we give a plot of the concurrence as a func- 
tion of temperature T for different anisotropic parame- 
ters. At zero temperature the concurrence is 1 since no 
matter what the sign of J is and what the values of 7 
are, the ground state is one of the Bell states, the max- 
imally entangled state. The concurrence monotonically 
decreases with the increase of temperature until it reaches 
the critical value of T and becomes zero. The numerical 
calculations also show that the critical temperature de- 
creases as the anisotropic parameter increases from to 
1. 
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FIG. 6. The concurrences versus temperature for different 
anisotropic parameters: 7 = (solid line), 7 = 0.6 (dashed 
line) and 7 = 0.8 (dotted line). The parameter J is set to 



VI. CONCLUSIONS 

In conclusion, we have presented some interesting re- 
sults in the simple XY model. First, we can use XY 
interaction to generate the 3-qubit and 4-qubit W entan- 
gled states. Second, we see that the time evolution of 
entanglement are periodic for 2, 3, 4 and 6 qubits, and 
there is no exact periodicity for large N. At some special 
points the states becomes disentangled. Finally we study 
the thermal entanglement within a two-qubit isotropic 
XY model with a magnetic field and an anisotropic XY 
model, and find that the thermal entanglement exists for 
both ferromagnetic and antiferromagnetic cases. Even in 
the simple model we see some evidence of the quantum 
phase transition. 

The entanglement is not completely determined by the 
partition function, i.e., by the usual quantum statistical 
physics. It is a good challenge to study the entanglement 
in multi-qubit quantum spin models. 
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